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C^ ' Abstract. We theoretically investigate pump-dump photoassociation of 

vj ' ultracold molecules with amplitude- and phase-modulated femtosecond laser 

ryy ' pulses. For this purpose a perturbative model for the light-matter interaction is 

^. ^ developed and combined with a genetic algorithm for adaptive feedback control of 

(-H , the laser pulse shapes. The model is applied to the formation of *^Rb2 molecules 

^ I in a magneto-optical trap. We find that optimized pulse shapes may maximize 

the formation of ground state molecules in a specific vibrational state at a pump- 
dump delay time for which unshaped pulses lead to a minimum of the formation 

fv^ ' rate. Compared to the maximum formation rate obtained for unshaped pulses 

^ , at the optimum pump-dump delay, the optimized pulses lead to a significant 

^^ ■ improvement of about 40% for the target level population. Since our model yields 

— ^ ' the spectral amplitudes and phases of the optimized pulses, the results are directly 

I ' applicable in pulse shaping experiments. 
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1. Introduction 

One of the challenges in contemporary atomic and molecular physics is the making of 
ultracold molecular samples ^2l,^.- Direct cooling techniques for molecules have not 
yet reached the ultracold regime (T < ImK). Instead, ultracold atoms are assembled 
to molecules using magnetic Feshbach resonances 4 or photoassociation |21 15] . In 
the latter case, typically a continuous-wave (cw) laser excites a pair of colliding atoms 
into a long-range molecular state. Subsequent spontaneous emission leads to molecules 
in their electronic ground state. In view of prospective applications, stable ultracold 
molecules, i. e. molecules in their absolute rotational, vibrational and electronic ground 
state are required. However, molecules formed with Feshbach resonances are created 
in the highest bound vibrational level T. Photoassociation with a cw laser followed 
by spontaneous emission forms molecules in a range of vibrational levels (2j. In 
this case, the radiative lifetime of the excited electronic state ultimately represents 
the rate-determining step in ground state molecule formation. When applied to a 
Bose-Einstein condensate (BEG), the incoherent nature of the spontaneous decay 
furthermore destroys the coherence of the BEG and hence the condensate itself |S]. 
This may be avoided by stimulated Raman transitions [H|. Recently, also light- 
driven Rabi oscillations between atoms and molecules have been observed |TI]] . Using 
cw photoassociation in combination with nanosecond pulsed Raman transitions, the 
formation of RbGs molecules in the v ~ level of the ground electronic state has been 
achieved {11 j . 

An alternative approach to ultracold molecule formation is the photoassociation 
with time-dependent laser fields of femtosecond to picosecond pulse durations. The 
conceptually simplest pulsed photoassociation scheme is based on a two-pulse pump- 
dump population transfer from the atomic continuum to a bound molecular level 
|12|. as shown in figure ^ The pump pulse excites a molecular wavepacket of the 
vibrational states within its bandwidth, which then propagates until it is dumped 
to the ground state by a second laser pulse. By modifying the delay time between 
the two pulses, the transition to one or several desired target vibrational levels in 
the ground state can be enhanced by quantum interference. Here, the production 
rate is not limited by the radiative lifetime time of the intermediate state, but by 
the much shorter vibrational period, by the Franck-Gondon factors and ultimately 
by the unitarity limit of the scattering rate jHl I13j . The particular advantage of this 
approach is that the pump-dump delay time and the pulse shapes yield additional 
control knobs to prepare specific product states, which is extensively explored in the 
field of adaptive feedback control with shaped femtosecond pulses ^^ E| . A number 
of theoretical studies have investigated the use of chirped laser pulses for excitation 
| 16l I17L [TH] and stabilization ^J [201 of ultracold molecules. Even larger control over 
the molecule formation process might be obtained by shaped laser pulses. Laser pulses 
which optimize the formation of ultracold molecules from an atomic BEG have been 
calculated using optimal control theory |2J| • In many cases the complexity of the pulse 
shapes resulting from optimal control calculations requires a careful transformation 
into experimentally achievable pulse shapes |22[ |23| . 

In this paper, adaptive feedback control of photoassociation with femtosecond 
laser pulses is investigated theoretically. For this purpose the properties of typical 
femtosecond oscillator pulses and a realistic pulse shaping device are combined with 
a closed loop genetic algorithm for pulse shape optimization. The resulting pulse 
shapes are obtained as amplitudes and phases in the frequency domain, similar to 
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Figure 1. The pump-dump scheme for the production of ultracold molecules 
with short laser pulses. A wavepacket is created in the excited state by a pump 
pulse (1) and subsequently propagates (2). A dump pulse transfers population to 
the target bound ground state level after a suitable time delay (3). 



the case in closed-loop adaptive feedback control experiments. The parameters of 
the calculations are chosen to correspond to a situation, which is typical for state- 
of-the-art experiments |24j . A two-channel model is employed, which simplifies the 
complex coupled-channel dynamics, but retains the important features of pump-dump 
photoassociation. The goal of this study is to qualitatively determine to which extent 
the formation process of ground state molecules can be controlled by realistic shaped 
laser pulses. 

A first experiment on photoassociation with femtosecond laser pulses has been 
carried out on hot mercury atoms j25| . In the ultracold regime, electronically excited 
Na2 molecules have been formed using picosecond pulses near-resonant to the sodium 
Dl transition, and probed with a subsequent laser pulse J5H]. Chirped nanosecond 
laser fields have been employed to enhance inelastic collision rates in an ultracold gas 
|27| . Also the effect of the chirp of femtosecond pulses on the dissociation of ultracold 
Rb2 molecules has been investigated |2H1- Recently, adaptive feedback control with 
shaped femtosecond pulses has been applied in our group to optimize the dissociation 
of ultracold Rb2 molecules from a magneto-optical trap i24j. We envision that the 
work presented in this paper will be helpful for the first experimental demonstration 
of ultracold ground state molecule formation using pulsed photoassociation, which has 
not been achieved to date. 

The paper is organized as follows: the model of two colliding atoms interacting 
with the laser fields is described in section 12.11 A perturbative treatment of the 
interaction with the field is employed in section 12.21 to obtain a model for the 
pulse shaping device. Its optimization with a genetic algorithm is outlined in 
section 12.31 The results are presented in section 13 The validity of the perturbative 
model is tested in section 13.11 by comparison to full, non-perturbative quantum 
calculations. In section 13.21 and section 13.31 the results for optimal laser pulses for 
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pulsed photoassociation, as obtained from adaptive feedback control are presented 
and analyzed. 

2. Model for the photoassociation process 

2.1. Interaction of two colliding atoms with a laser field 

A colliding atom pair subject to a laser pulse is considered. Two electronic states 
coupled by the radiation field are taken into account with the assumption of the 
rotating wave approximation. The center-of-mass motion of the two atoms is assumed 
to be decoupled from the internuclear dynamics. The time-dependent Schrodinger 
equation then reads 

^^^d_r yjg{R,t) \_f,f Vg(i?,t) 
with 



'^ dt( MR,t) ) ^[ MR,t) ' ^^^ 



\ hn*{t) f + Ve{R)-A J ' ^ ' 

The wavefunction consists of one component for the ground state, '4>g{R,t), and one 
component for the excited state, ipeiRjt). The diagonal elements of the Hamiltonian 
are given by the kinetic energy operator T and the potential energy curves for the 
respective states Vg/f.{R) which depend on the internuclear separation R. The zero 
of energy is chosen to correspond to the ground state dissociation limit. Within 
the rotating wave approximation, Ve is shifted by the detuning A with respect to 
Vg where A = h{uige — ojl)- ^ge denotes the atomic transition frequency and uj^ 
the carrier frequency of the laser. The off-diagonal elements of the Hamiltonian are 
given by the complex time-dependent Rabi frequency ri(t), which is related to the 
electric field S{t) by Q,{t) = ij,£{t)/fi. Here, /i denotes the dipole matrix element 
of the two molecular electronic states under consideration. In the following, the in- 
dependence of this matrix element is neglected and its asymptotic value is employed, 
fi = 2.537 • 10~^ Cm for the Rb Dl line. The laser pulse is assumed to be derived from 
a Gaussian transform-limited pulse sent through a pulse shaping device as detailed in 
section lOl 

The two channels under consideration for the calculations are the singlet ground 
state X^Eg(5s -I- 5s) and the excited state 0+(5s -|- 5pi/2) of *^Rb2. The wavelength 
of the asymptotic (atomic) transition corresponds to 794.9 nm. The potential energy 
curves are derived from ab initio data |29| at short distances and long-range dispersion 
potentials {—C^/R^) — Cq/R^ — Cg/R^ for distances larger than the LeRoy-Bernstein 
radius. The coefficients for the 5s -I- 5s asymptote are found in jSO], while the 
coefficients for the 5s -I- 5p asymptote are taken from jSJ- The Hund's case (c) 
representation of the potential energy curves is obtained by diagonalizing the Hund's 
case (a) potentials coupled by spin-orbit interaction [32|. The spin-orbit coupling 
matrix element is set to the atomic value. The i?-dependence and the resonant 
character of the spin-orbit coupling are neglected. The latter leads to strong 
perturbations in the vibrational spectrum |83| providing an efficient stabilization 
mechanism for pump-dump photoassociation |2()| . However, its treatment would 
require a three-channel model, and the focus of this study is on the qualitative 
understanding of possible control over the wavepacket dynamics. 
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The Hamiltonian (jSJ is represented on a grid using the mapped Fourier grid 
method OOS]- A grid of 1024 points extending up to 27000 Bohr radii is chosen. 
This size is necessary to obtain a reliable representation of continuum scattering states 
| 18| . The binding energies i^j? and wavefunctions (pv (R) of the vibrational levels 
are computed by diagonalizing the Hamiltonian (O with n{t) = 0. Franck-Condon 
factors are then obtained as scalar products, |(?;|w)p = / dRlipKR)* ipf^{R)\'^ . To test 
the validity of the perturbative model which is explained in the next section, the time- 
dependent Schrodinger equation ^ is solved with the Chebyshev propagator method 

The initial state of the colliding atom pair is chosen such as to describe an 
ultracold atomic gas in a magneto-optical trap. Due to the low temperatures of about 
100 /zK, only a narrow band of continuum states above the ground state potential 
asymptote is occupied. The considerations in this study focus on a single initial 
continuum state with coUisional angular momentum I — 0, i.e. an s-wave scattering 
state. This assumption is justified since the continuum wavefunctions corresponding 
to temperatures T < 500 /iK possess an identical nodal structure at short and 
intermediate internuclear distances. A thermal average over all possible initial states 
is only required for absolute molecule formation rates, which is beyond the scope of 
the present study. 

2.2. Model of the photoassociation process with shaped laser pulses 

In the following a perturbative model of pump-dump photoassociation with arbitrarily 
phase and amplitude modulated laser pulses is presented. Such shaped laser pulses are 
prepared experimentally by passing the spectrally dispersed pulses through a liquid 
crystal modulator array. In each pixel of such a pulse shaper the phase and the 
amplitude of the corresponding frequency component of the pulse can be computer 
controlled. 

The photoassociation pulse creates a time-dependent wavepacket in the excited 
state. After the pulse is off, the excited state wavefunction ipeiR'jt) is given by a 
superposition of vibrational levels in the excited state. If the interaction can be 
treated perturbatively, i.e. the Rabi angle, — J il(i) dt, is small compared to tt, 
the pulse excites the superposition of vibrational levels given by |87j 

V'e(i?;0«I]a. ^('-'«)('^l^)</'.(^)e-'''"*^^ (3) 

V 

where {v\E) is the overlap integral of the initial continuum state \E) and the excited 
state vibrational level w; its square corresponds to the free-bound Franck-Condon 
factors. The frequency spectrum of the electric field of the laser pulse, S{uj), is 
evaluated at the transition frequency uj^ to the vibrational level v. The a^ coefficients 
are introduced to represent the complex modulation coefficients of a pulse shaper. 
This description assumes that the frequency components for each contributing v are 
individually resolved by the pulse shaper. As one can experimentally only diminish 
spectral components, one obtains the constraint < |at,| < 1. 

Equation Q is only valid, if the pulse does not contain significant amplitude at 
frequencies close to resonance with the atomic transition. This is due to the strongly 
increasing transition strength close to the atomic resonance, which leads to a break- 
down of the perturbative treatment for typical laser intensities. We estimate that the 
perturbative model is valid if the detuning of the laser pulse spectrum with respect 
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to the atomic transition is large enough so that minimal excitation of the atomic 
transition occurs. After a delay time At, a second or dump pulse £'{uj) is applied. 
This pulse transfers population from the excited state to the ground state; it is assumed 
to be not shaped by a pulse shaper. The resulting superposition of bound ground state 
levels reads 



-i-_Bf„(t-At) -i-B^At 



{E\v){v\w) (^^(i?)e-«^™^*-^*Je-^^"^\ (4) 

Note that this description is only valid for non-overlapping pulses, because transient 
dynamics that occur during the excitation process are not included. Therefore only 
delay times At that are large compared to the duration of the pump and dump pulses 
are chosen. In Q, u runs over all bound vibrational levels of the excited electronic 
state and w over all bound levels of the ground state. The {v\w) are overlap integrals 
between ground and excited state bound levels. The time delay between the two pulses 
enters the expression through the phase factors ^-'^^l^^lf^ which corresponds to the 
propagation of the wavepacket in the excited state. The population of an individual 
bound ground state level w is now given by 

v,v' 

■ {E\v){v\w){w\v'){v'\E)e-^'^^i'-^^^'^^K (5) 

This double sum over the excited state vibrational levels v and v' can be concisely 
written as a complex quadratic form in terms of the pulse modulation coefficients ay , 

|6^P = ^a:, Mr'(At)a.. (6) 

The Franck-Condon matrix M™ (At) contains the properties of the system (overlap 
integrals and binding energies) and of the frequency spectrum of the two laser pulses 
and their pulse delay. 

It can be inferred from that the choice of modulating the pump pulse and not 
the dump pulse is in fact arbitrary, since the same structure of is obtained if the 
dump pulse or both pulses are shaped. Equation (jB)) serves two purposes: First, for 
a given target level w, it can be used to obtain the optimal detuning of the pulses by 
evaluating |6u,p for different delay times At. Second, it allows for the implementation 
of an optimization procedure as explained in the following section. 

2.3. Genetic algorithm for the pulse optimization 

In the following, the photoassociation pulse shall be optimized while the dump pulse is 
assumed to be transform-limited. An optimization algorithm is obtained by combining 
the analysis of the target state © with a spectral decomposition of the pulse. 
The photoassociation pulse which yields the desired excited state wavepacket can 
in principle be spectrally composed from a set of weighted delta functions at the 
transition frequencies corresponding to the vibrational states v |37| . This spectrum of 
a few sharp structures leads to very long pulses in time domain, which is unfavorable 
for our perturbative model of a pump dump scheme as explained above. The 
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photoassociation pulse is instead constructed from a coherent superposition of a larger 
number of frequency components, 

hn{t) = ^iY. c„f K)e^(""-"^)* , (7) 

n=l 

which leads to a shorter pulse in the time domain. The frequency spectrum of the 
pump and dump laser pulses before modifications by the pulse shaper is assumed to be 
a Fourier limited Gaussian pulse £{ujn) oc exp [— (w„ — lul)'^ /2T^] with a root-mean- 
square width r. Each a;„ corresponds to a pixel in the pulse shaper mask, lol is again 
the carrier frequency of the laser pulse. The pixel number Np is set to 640 and the 
frequency spacing from pixel to pixel is set to 1.5 cm~^, which are numbers used in 
experimentally available pulse shaping devices. This leads to a pulse shaper resolution 
that is smaller than the vibrational level spacing for almost all the excited vibrational 
levels. The pulse shaper coefficients c„ are set equal to the desired superposition 
coefficients a„ from ©, where v is the excited state level for which hE^ is closest to 
LOn- The range of vibrational levels which contribute to the population in the excited 
state is identified by the non- vanishing matrix elements of Af ™ . In the present case, 
this range consists of 150 levels ranging from v — 170 (Tvib = 0.77 ps) to w = 320 
(Tvib = 80.0 ps). The discreteness of the spectrum ((TJ still leads to unwanted replica 
effects in the time domain. The pulse is therefore multiplied with a box function of 
the form exp [— (t/Tcutoff)™] with m = 8 and Tcutoff' = 3.33 ps; t is the time relative to 
the maximum of the pulse envelope. This provides a smooth switching on and off of 
the pulse. 

The optimization is carried out with respect to the a„ in (jSJ, corresponding to 
the pixels of a pulse shaping device. The pump-dump delay time At represents a 
fixed parameter for the optimization process. It turns out that for the amplitude 
optimization it is sufficient to set a„ to either or 1 . For the phase 256 possible values 
are considered, leading to an overall 9-bit optimization of the a^ . The quality of the 
pulses is determined by the fitness factor Q which is chosen as 



Here, (i^) is a measure of the pulse duration, 

^ ^~ jn{t)dt ■ ^ ' 

It is included in the definition of the fitness factor in order to keep the pulses as short as 
possible. This is not a physical constraint but a requirement of the perturbative model, 
because the duration of the pulses is immensely increased by the modulation: Already 
a linear frequency chirp stretches the pulse in time. The employed perturbative model, 
however, requires non-overlapping pump and dump pulses, as explained above. 

The optimization of the formation of molecules in a specific ground state level 
w is carried out by means of an evolutionary algorithm. This corresponds directly 
to a close-loop experiment. For the calculations a very basic evolutionary algorithm, 
SimpleGA from the GALib package |3H1, is employed. Populations of 20 individuals 
are used, and the number of generations is set to 1000. This is sufficient to obtain 
convergent results. The single-bit mutation and single-bit crossover probabilities are 
chosen as 0.005 and 0.8, respectively. 
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Figure 2. Dependency of the excitation probability on the peak intensity of the 
pump pulse for a fixed pulse width of 100 fs FWHM. The linear dependence on 
intensity below 2 MW/cni'^ shows the validity of a perturbativc treatment. 



3. Results 



3.1. Validity of the pulse shaper model 

The pulse shaper model of © is derived for the case of a small perturbation of the 
atom pair by the laser field. Higher order processes, in particular Rabi cycling during 
the pump or dump pulses, are neglected. By carrying out calculations for a single 
unshaped pump pulse with varying the peak intensity of the pulse, the regime of 
validity of the perturbativc approximation is analyzed. 

In figure |2] the probability to excite an atom pair to the excited state is shown 
as a function of the pulse's peak intensity. These results are obtained by solving the 
time-dependent Schrodinger equation where the field is treated non-perturbatively. 
Figure [2 clearly shows a linear dependence of the excitation probability for small 
intensities < 4 • 10® W/cm^. For these low intensities the effect of the laser field can 
be treated in first order perturbation theory. For larger pulse energies the excitation 
starts to saturate, which indicates the transition to a non-perturbative regime. This 
is attributed to Rabi cycling and strong off-resonant excitations. 

In order to check the validity of © in the perturbativc regime, a model study is 
performed. The excitation of a vibrational wavepacket in the excited state, tpe [R] t) 
is calculated with ||2J) for a single Fourier-limited Gaussian pump pulse. The detuning 
of the pulse is set to A = 251 cm"-'^ and the temporal intensity width to 100 fs 
FWHM. These values match the experimental conditions of Ref. [23 ■ In figure|21 the 
projections of the excited state wavepacket onto the excited vibrational eigenfunctions 
(fil {R) are shown as a function of the detuning, E^ (solid line) . This is compared in 
figure 121 to the non-perturbative calculation for the same pulse parameters (squares). 
Good agreement of the projections is observed in the range of resonant levels. 

As described above, the perturbativc pump-dump model is only valid for non- 
overlapping pulses, because transient dynamics occurring during the pulsed excitation 
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Figure 3. Comparison between tlie predicted excited state level populations 
after the pump pulse from the model (solid line) and the results of the quantum 
dynamical calculation (squares). 



are not included. Therefore, in a second step it is checked if © holds for two-pulse 
pump-dump processes. Assuming two pulses with width and detuning, as given above, 
the pump-dump process is calculated with 101 for varying delay time. This yields the 
ground state wavepacket which is then projected onto a single target vibrational level. 
The level w = 119, bound by Ef^g =6.86 cm~^, is selected, because it features a large 
Franck-Condon factor from the excited state vibrational levels that are populated 
at the given detuning of the laser pulse (see section 13.21 and figure O below) . The 
population in this target level is shown in figure ^ (solid line) as a function of the 
pump-dump delay time. It is compared to the results obtained by solving the time- 
dependent Schrodinger equation (squares). The agreement between the perturbative 
model and the full quantum dynamics with a non-perturbative treatment of the field 
is convincing. 

Note that in the present case, the optimal pump-dump delay is given by the full 
excited state wavepacket round-trip time Tvib and integer multiples thereof. According 
to the Franck-Condon principle, this indicates that the dump transition occurs at the 
outer turning point of the excited state wavepacket, as indicated in figure ^ This is 
not necessarily the case, since it depends strongly on the molecular potential curves 
and the laser detuning. The dump transition will occur at half-integer multiples of 
Tvib, corresponding to the inner turning point if the potential facilitates a deceleration 
mechanism at short distances. This is the case for a soft repulsive wall such as in the 
outer well of the 0~(p3/2) state ^^Ij or for resonant coupling |20| . 



3.2. Optimal detuning 

In order to find a window of detunings suitable to populate target levels in the ground 
state, © is evaluated for unshaped Fourier-limited Gaussian pulses of 100 fs FWHM 
intensity width. For a range of delay times encompassing several round-trip times, the 
Franck-Condon matrix M™ (At) is calculated for each ground state target level w. 
The optimal detuning is inferred from the maximunr of the main diagonal of Af ™ ( Ai) 
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Figure 4. Target state population |6™p for w = 119 from J^J (solid line) 
compared to the results of the quantum dynamical calculation (squares). The 
peaks are closely connected to the wavepacket dynamics in the excited state 
potential well. The temporal smearing out of peak heights is due to anharmonic 
dispersion of the wavepacket. 
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Figure 5. Optimal detuning (black ■, left axis) for pump-dump photoassociation 
with two identical unshaped femtosecond laser pulses, estimated from j^J for a 
range of target vibrational levels w in the ground state. The maximum attainable 
population of these levels (red • ) is shown on the right axis. 
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and the corresponding maximum target state population is calculated. 

FigureElshows the result for both the optimal detuning and corresponding target 
state population as a function of the vibrational level w of the ground electronic state. 
A strong dependence of the maximum attainable population on the vibrational level is 
observed spanning eight orders of magnitude. The largest target populations are found 
for small laser detunings and correspond to high target vibrational levels. This reflects 
the large free-bound Franck-Condon factors for highly excited vibrational levels. It 
also shows that the bound-bound Franck-Condon factors for the dump step favor high 
target levels for small detunings. 

3.3. Adaptive feedback control 

With the perturbative model, optimization calculations are carried out for pump- 
dump photoassociation. As in the previous section, the pump and dump pulses are 
assumed to be Fourier-limited Gaussian pulses with 100 fs FWHM intensity width 
and A = 251 cm~^ detuning. The goal for the optimization is to find the pump pulse 
that maximizes the population in the w = 119 ground state vibrational level, which 
is a suitable level for the chosen detuning due to its large Franck-Condon overlap (see 
section m|l . The spectral components of the pump pulse are optimized with respect to 
their phase and amplitude. The pump-dump delay time, which represents a parameter 
in these calculations, is set either to 15.75 ps or to 14.0 ps: As shown in figure 01 for 
At = 15.75 ps a maximum is observed in the target level population of the 'natural' 
dynamics, i.e. the dynamics obtained with an unshaped pump pulse. On the other 
hand. At = 14.0 ps represents a delay for which the natural dynamics leads to almost 
no target level population. 

The resulting target state population with optimized phases and amplitudes is 
shown as a function of the pump-dump delay time in figure (solid lines) . Note 
that the target level population on the vertical axis of both graphs in figure El and of 
the graph in figure^ is given in the same units and can thus directly be compered. 
In figure it is clearly seen that for both optimization calculations the maximum 
population is indeed obtained at the delay time for which the pump pulse is optimized. 
For At — 14.00 ps, the dynamics is modulated such that an efficient dump transition 
may occur at delay times when almost no target level population is obtained with 
an unshaped pump pulse, see the left arrow in figure ^ In this case, the target 
level population for the optimized pulse amounts to 1.6 arb. units., whereas less 
than 0.1 arb. units are observed for the same pump-dump delay with an unshaped 
pulse, corresponding to an increase by increased by more than a factor of 10. This 
illustrates that the pulse shaper can be used to a great extent to control the wavepacket 
dynamics. The optimal pulse shapes also lead to a higher maximum of the target level 
population as compared to the maximum population achieved for the optimal pump- 
dump delay time with unshaped pulses: for the optimization at At — 14.0 ps, the 
target population is enhanced by about 40%. For 15.75 ps a smaller enhancement of 
about 10% is observed. The two different levels of optimization are attributed to the 
different integrated energies of the obtained optimal pulse, as detailed below. 

The results of the optimization are compared with the results of the time- 
dependent Schrodinger equation using the optimal pulses, obtained above, as pump 
pulses. The time delay between the pump and dump pulses is varied to investigate 
whether the optimal result is indeed obtained for the intended delay and how in general 
the pulse shapes affect the excited state wavepacket dynamics. As shown in figure 
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Figure 6. Relative target state {w = 119) population versus pump-dump delay 
for optimized pump pulses. The target level population on the vertical axis of 
the graphs is given in the same units as in figure |1] and can thus directly be 
compered. The optimizations are carried out with a delay parameter of 14.00 ps 
(upper trace) and 15.75 ps (lower trace), corresponding to a minimum and a 
maximum in the 'natural' dynamics after an unshaped pump pulse, cf. figure El 
The pulse parameters are the same as in figure |2l The prediction of the model 
151 (solid line) is compared to the solution of the time-dependent Schrodinger 
equation (squares). Note the deviations occurring at small delay times which are 
due to transient excitations. 



these results agree well with the model, except for very short delay times where the 
pulses still overlap. 

The mask patterns of the pulse shaper for amplitude ( ) and phase ( ) 

are shown for the two preset delay times in the top panels (a) and (b) of figure [3 
The spectrum contains frequencies for which the excited state levels v = 204 to 
V = 249 are resonantly excited from the initial state. For At = 14.0 ps, figure 
(a), almost the full spectrum is retained, whereas for At = 15.75 ps, figure[3(b), some 
frequency components are removed from the pulse spectrum. This leads to a smaller 
integrated energy of the optimized pulses for the At = 15.75 ps optimum. Therefore a 
smaller target level population is observed on maximum (about 1.3 arb. units) for the 
15.75 ps calculation as compared to the 14. Ops calculation (1.6 arb. units). Chirp-like 
structures in the spectral phase emerge in both cases. In other runs of the optimization 
algorithm, such structures have been obtained as well. 
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Figure 7. Upper panels: Amplitude (red ) and phase (black ) mask 

patterns of the pulse shaper to generate the optimized pulses for At = 14.00 ps (a) 
and At = 15.75 ps (b). The amplitudes are multiplied by the spectral envelope of 
the pulse. Note the chirp-like structures in the phases of the two pulses. Lower 
panels: Temporal amplitude of the optimized pulses for the 14.00ps case (c) and 
the 15.75ps case (d). 



The temporal profiles of the optimized pulses are shown in the bottom panels (c) 
and (d) of figure The pulses are stretched in time compared to the time duration 
of the unshaped pulses of 100 fs FWHM. For At — 14.0 ps, figure [7|(c), three sub- 
pulses are observed which are separated by about 2ps. This has to be compared to 
the vibrational times of the excited state levels which range from 2.7 ps for v = 204 
to 10.0 ps for V = 249 and 4ps for the central level v = 219 resonant with the carrier 
frequency. In the Wigner representation of the optimal field (data not shown), an 
interference pattern between the three sub-pulses is observed indicating the intricate 
phase relation between the sub-pulses. The 'holes' in the spectrum observed for 
At = 15.75 ps, as shown in figure[71(b), correspond mostly to transition frequencies to 
levels with long vibrational periods (for example, v = 227 to w = 233 with Tvib = 5.1 ps 
to 6.3 ps). 

Certain wavepacket components are filtered out to facilitate a shaping of the 
excited state wavepacket which is favorable for the dump step. This is illustrated for 
At — 15.75 ps in figure |H1 The excited state wavepacket at t = 15.75 ps after the 
unshaped, figure |H1 (a), and optimized, figure |Hl(c), pump pulse is compared to the 
vibrational wavefunction of the ground state target level, figure|Hl(b). Optimal dump 
conditions are obtained when the wavepacket is focused in the region where the target 
state wavefunction possesses its outermost maximum. The optimization leads to a 
removal of parts of the wavepacket in the inner region which would cause destructive 
interference in the dump process. 
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Figure 8. Excited state wavepackets in position space at t = 15.75 ps after the 
pump pulse maximum. The upper panel (a) shows the wavepacket obtained after 
an unshaped pulse, while the lower (c) panel displays the wavepacket resulting 
from excitation with a pulse optimized for At = 15.75 ps. To estimate the dump 
efficiency, these wavepackets have to be compared to the stationary ground state 
wavefunction of the target level shown in the middle panel (b) . Note the reduced 
probability density at small internuclear distances of the wavepacket after the 
optimized pulse as compared to the unshaped case. 



4. Conclusions and outlook 



We have investigated to which extent the formation of ultracold molecules by pulsed 
photoassociation can be controlled using adaptive feedback control and shaped 
femtosecond laser pulses. For this purpose we have set up a perturbative pulse 
shaper model and combined it with a genetic algorithm. This concept allows for 
a numerically efficient optimization with realistic conditions for the femtosecond 
laser pulse and the pulse shaper. Using the model for a pump-dump scheme, the 
population in a specific target level in the ground electronic state is optimized. 
A considerable enhancement factor of about 40% is achieved compared to highest 
population achieved with unshaped pulses, while keeping the pulse duration as short 
as possible. The maximum pump-dump enhancement in the perturbative, low-energy 
regime, is determined by the Franck-Condon factors of the two electronic transitions. 
Significantly larger enhancements may only be expected for high-intensity pulses 
beyond the perturbative regime, where Rabi cycling and dynamic Stark shifts occur. 
The adaptive feedback optimization has shown that for different time delays 
between the pump and dump pulses, a maximum in the ground state molecule 
formation rate of similar height is achieved. This implies that the pump-dump delay 
represents an uncritical parameter as long as it is shorter than the spontaneous decay 
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time, in the sense that the optimization algorithm may find an optimal solution with 
a comparable target level population for arbitrary delay times. The analysis of the 
excited state wavepacket dynamics reveals that the shaped pulse induces a spatial 
focusing of the wavepacket in the Franck-Condon region of the dump pulse at the 
specified pump-dump delay. The pulse shaper can therefore greatly alter the dynamics 
of the wavepacket to optimize ground state molecule formation. 

The model described in this work is advantageous compared to standard 
theoretical approaches such as optimal control theory [21] in that it closely resembles 
the situation in experiments. Optimal control calculations generally converge faster 
than genetic algorithm optimizations and extremely efficient pulses can be obtained 
PO). However, these pulses are usually too complex to be directly implemented in an 
experiment J22II28I . In this work, the optimized pulses are obtained in the frequency 
domain and are given in terms of the transmission coefficients and phase shifts of 
a pulse shaper. The optimization results can therefore directly by be employed in 
experiments. 

Future improvements of the presented calculations may follow several lines. To 
obtain absolute rates for molecule formation, a thermal average needs to be performed. 
This will be important to find pulse shapes for which association dominates over 
photodissociation of ultracold molecules |24l I28j . For the collision energies in a 
magneto-optical trap, such an average requires inclusion of partial waves with I > 0. 
Second, all relevant electronic states need to be incorporated in the model. In 
particular, resonant spin-orbit coupling is known to lead to perturbations in the 
excited state vibrational spectrum |H8I I41| which might be important for ground state 
molecule formation 20 . Furthermore, it will be interesting to investigate the influence 
of the laser field polarization on ground state molecule formation, in particular to 
find shaped pulses which excite specific molecular symmetries. Finally, higher order 
terms may be included in the perturbative derivation of the pulse shaper model to 
allow for investigation of the nonlinear excitation regime. Here, higher enhancement 
factors beyond the Franck-Condon regime may be expected. Based on the findings 
in this work, dedicated experiments on pulsed photoassociation of ultracold rubidium 
molecules are in progress. 
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